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Abstract 

Approximating the sum over all gauge field configurations in the QCD partition func- 
tion by a liquid of instantons, we calculate the spectrum of the Dirac operator for two and 
three colors and for 0, 1 and 2 flavors. We find a remarkable difference in the spectrum 
near zero virtuality between 2 and 3 colors, which can be explained in terms of chiral ran- 
dom matrix theory. For two colors the Dirac operator is real, and the appropriate random 
matrix ensemble has real matrix elements. For three colors the Dirac operator is complex, 
and the spectrum can be described by a random matrix ensemble with complex matrix 
elements. These results provide further evidence that the spectrum near zero virtuality 
is universal and is completely determined by symmetries. 
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1. Introduction 



It is widely believed that at some critical temperature the phase of QCD changes 
from a chirally broken phase to a chirally symmetric phase. The order parameter is the 
chiral condensate, which, according to the Banks-Casher formula [ffj, is proportional to 
the thermodynamic limit of the spectral density of the Dirac operator at zero virtuality. 
The study of this order parameter is the central objective of this paper. In particular, we 
will analyze the spectrum near zero virtuality via the microscopic limit which is obtained 
in the thermodynamic limit V4 — > 00 by rescaling the eigenvalues by V4. 

A natural way to perform these studies would be via lattice QCD calculations. How- 
ever, the order parameter is defined in the chiral limit which makes such calculations with 
dynamical fermions infeasible. Instead, we approximate the integral over all gauge field 
configurations by a liquid of instantons and anti-instantons. In this approach the chiral 
limit does not give rise to technical problems. 

Our main goal is to substantiate the claim || that spectral correlations on the scale of 
one average level spacing are universal. It is well known from the study of systems that are 
classically chaotic ||, ||J] H, || [0 that correlations between eigenvalues measured in units 
of the average level spacing are universal in the sense that they are solely determined 
by the symmetries of the system. Therefore they can as well be described by random 
matrix theory with only the symmetries as input. It is our conjecture that the same is 
true for correlations of the eigenvalues of the Dirac operator measured over the ensemble 
of all gauge field configurations. In QCD, though, we have an important additional 
symmetry: the chiral symmetry. As a consequence all nonzero eigenvalues occur in pairs 
±A. Obviously, this symmetry has to be implemented in the random matrix model. 
What is important is that for these types of random matrix models, which we call chiral 
random matrix models, the point A = is special: the spectrum has a hard edge at A = 
II |!| . Therefore, correlations on the scale of one average level spacing can be defined 
with respect to this point. The microscopic spectral density is exactly such correlator 
and is expected to be universal. There are other arguments in favor of universality. Most 
notably, it was recently argued by Leutwyler and Smilga |i| that the low-energy effective 



partition function of QCD which is solely based on symmetries leads to an infinite family 



of sum rules for the inverse powers of the eigenvalues of the Dirac operator. Since all 
sum rules can be derived from the microscopic correlators, this strongly suggest that the 
spectrum at this scale is determined by symmetries as well. 

For the classical random matrix ensembles we have three universality classes: the 
GOE, the GUE and the GSE, depending on whether the matrix elements are real, com- 
plex or quaternion real, respectively [|TTJ. The same is true for the chiral random matrix 
ensembles. Recently, we succeeded the classify the SU(N C ) gauge theories accordingly 



J T2] . For N c = 2 and fundamental fermions all matrix elements of the Dirac operator can 
be made real, for N c > 3 gauge groups with fundamental fermions the Dirac operator is 
complex, whereas for adjoint fermions, the Dirac operator can be regrouped into quater- 
nions. The corresponding random matrix ensembles will be called the chiral orthogonal 
ensemble (chGOE), the chiral unitary ensemble (chGUE) and the chiral symplectic en- 
semble (chGSE), respectively. Each of the three cases can be identified with one of the 



three possible scenarios of chiral symmetry breaking [14]. Up to now, the microscopic 



spectral density has been obtained for the chGOE |15j and the chGUE |]13 



As a final note in this introduction, we want to point out that the microscopic fluc- 
tuations become universal only after they are separated from all other variations in the 
spectrum. This procedure, called unfolding |§, will be applied to all our numerical results. 
It can only be achieved if a separation of scales takes place. We want to stress that this 
implicit assumption is a nontrivial point, that cannot be taken for granted. 

In this paper we calculate the spectrum of the Dirac operator for SU{2) and 577(3) 
color gauge groups and compare the results with the universal functions obtained from 
random matrix theory. The microscopic spectral density is defined in section 2. The main 
ingredient of the instanton liquid model are discussed in section 3. The expressions for 
the microscopic spectral density are given in section 4. Numerical results are presented 
in section 5 and concluding remarks are made in section 6. 



2. The microscopic spectral density 

For Nf flavors with masses rrif (rrif — > 0) the QCD partition function in the sector 



3 



with v zero modes is defined by 

2? cD H = <n nw+^/Kkw- (i) 

/=1 A n >0 

where the average is over gauge field configurations with v fermionic zero modes weighted 
by the gauge field action S„(A). The product is over all eigenvalues of the Dirac operator. 
The spectral density of the QCD Dirac operator for a given field configuration is defined 
by 

p(X)=J2S(X-X n ). (2) 

n 

The average spectral density (p(A)) is defined with respect to the partition function (1). 
According to the Banks-Casher formula jl| the chiral order parameter can be expressed 
as 

<^> = lim lp(0), (3) 
V4-+00 1/4 

where it is understood that the limit V4 — > 00 is taken before the chiral limit m — > 0. 
Its existence implies that the eigenvalues near zero scale as I/V4. Therefore, the order 
parameter can be analyzed in much greater detail via the microscopic limit of the spectral 
density defined as 

Ps(z) = lim I )\ . (4) 

y 4 ^oo k 4 \ 1/4/1/ 

Instead of a single number, our new order parameter is a function that not only gives the 
thermodynamic limit of (ipip)-, but also provides information on how the thermodynamic 
limit is approached. 



3. Instanton liquid model 

The sum over all gauge field configurations in (1) has only been performed exactly 
in lattice gauge theory. However, in calculations with dynamical fermions, the quark 
masses cannot be taken equal to zero, which makes it difficult to study the chiral limit of 
the microscopic spectral density. This problem is not present in our approach, where we 
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approximate the sum over all gauge field configuration in (1) by a sum over instantons 
and anti-instantons. 

In this approximation, the non-zero mode quantum fluctuations about the classical 
configurations are taken into account to one-loop order under the assumption that con- 
tributions from different instantons factorize. This is justified when the ensemble of 

is 



instantons is sufficiently dilute, which according to phenomenological arguments |L6| 
indeed the case. The integral over the zero modes, i.e. the size, the position and the 
orientation of each of the pseudoparticles, is done exactly via a Monte Carlo simulation. 
Also the integral over the fermionic zero modes is performed exactly. The latter amounts 
to evaluating the fermion determinant in the space of fermionic zero modes under the as- 
sumption that it factorizes from the contribution of the non-zero modes. Again, this can 
be justified on the basis that the instanton ensemble is sufficiently dilute. The resulting 
Dirac operator is given by the (Nj + Na) x (Nj + Na) matrix 



irrif T 
T' irrif 



(5) 



where the overlap matrix elements are given by 

T IA = [ d^ 1 (x)W^(x) = 1 Tr (t+R^U^Ua) F(X), (6) 
J APi Pa) 2 

Here, Rj, Ra are the positions of the instanton and the anti- instanton, pj, pa, their sizes 
and Uj and Ua their orientations in color space. The vector R IA is a unit vector in the 
direction (Rj — Ra)^, and we have introduced the notation r+ = (r, i)^ with f are the 
Pauli spin matrices. The scalar function F(X) depends on the specific ansatz for the gauge 
field configuration. We use the so called streamline field configurations |T7], [HJ, which 
are optimal in the sense that the first derivative with respect to directions perpendicular 
to the collective coordinates vanishes. In this case A is given by a conformal invariant 
combination of Ria, Pi and pa- Asymptotically, we have 

F(A)~4%^ (7) 

n IA 

At this moment we want to point out that for SU(2) all matrix elements of T are real, 
whereas for SU(3) and all larger gauge groups they are complex. As we will see below, this 
distinction has important consequences for the spectral correlations near zero virtuality. 



Summarizing, the QCD partition function is approximated by 

N N f 

Y[ dUidzidpifi^piA) Y[ 

i=l f i<j 



N N f 

J J] dUidZidp&ipiA) IJ det(TtT + mj) exp(-/3( P A) £ Sff), (8) 



where p(pA) is size distribution as induced by the quantum fluctuations (A is the QCD 
scale parameter). The total density of instantons N/V4 is kept fixed (N = Nj + Na), 
and we always work in the sector of zero total topological charge. For the interaction 
of the instantons we use the streamline action supplemented by a hard core. The latter 
has been introduced in order to be able to simulate a stable liquid of instantons. In this 
way, instantons are simulated as they occur in the full QCD partition function as given 
by lattice QCD simulations []19||. For a more detailed discussion of the instanton liquid 



model described above we refer to refs. 20, 21, 22 



4. Random matrix theory 

Our hypothesis is that the fluctuations of the eigenvalues of the Dirac operator near 
zero virtuality are universal and can be described by chiral random matrix theory. Using 
the instanton partition function as an inspiration p3|, |24|, M it is clear that the appropriate 



random matrix model in the sector with v fermionic zero modes is defined by [I2| 



J /3,v 



m) 



/«rw>3d*($ £), (9) 



where T has the symmetries of the corresponding Dirac operator and the masses are in 
the chiral limit (rrif — > 0). As discussed above, depending on the universality class the 
matrix T is real ((3 = 1, chGOE), complex (f3 = 2, chGUE) or quaternion real (j3 = 4, 
chGSE). In the latter case the square root of the fermion determinant appears in (9). 
The matrix T is a rectangular n x m matrix with \n — m\ = v (for definiteness we take 
m > n), so that the 'Dirac operator' in (9) has exactly v zero modes. The function P(T) 
is chosen gaussian: 

^) = exp(-^£A^. (10) 



In this normalization, the chiral condensate (as given by (3)) equals S in each of the three 
random matrix ensembles. The total number of modes is N = m + n. The latter quantity 
is identified with the volume of space time. 

The spectral density can be obtained by transforming to new integration variables in 
which T is diagonal. For (3 = 2 the corresponding Jacobian can be rewritten in terms 
of generalized Laguerre polynomials. After some manipulations, which are well-known 
from random matrix theory || , one obtains a formula for the spectral density in terms of 
these polynomials. From their asymptotic properties it follows that the microscopic limit 
is given by 

Y?z 

Ps{z) = —(JN f+ v(^z) ~ Jjv /+ i/+i(Sz) Jjv>+v-i(IU0). (11) 

In a more modern language, this limit is a double scaling limit, which exists as a conse- 
quence of the fact that the Laguerre ensemble has a hard edge at zero || |J . 

The situation for the chGOE is much more complicated. Because of the structure of 
the Jacobian, the standard orthogonal polynomial method no longer works. However, 
in this case it is possible to express the spectral density in terms of skew-orthogonal 
polynomials p5f . Remarkably, using a number of tricks from f26[ it is also possible in this 



case to obtain an analytical result for the microscopic spectral density |0|. The result 
can be expressed as an integral over the Bessel kernel 



Ps{z) = -J 2a+ i(zE) + - , ; i 

4 2 Jo \w aw z dz / 

wJ 2 a{zT)J 2 a~l{wT 1 ) - zJ 2 a~l(zT l )J 2 a(wT,) 



(zw) 2a (z 2 — w 2 ) 



(12) 



The parameter a is defined by 



a = N /~\ + ? ( 13 ) 

Recently, the Bessel kernel got a great deal of attention in the literature. For details and 
further references we refer to M. PI . 



5. Numerical results 
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All results have been obtained for an ensemble of 32 instantons and 32 anti-instantons 
confined to a box of 2.37 3 x 4.74A^ D . For N f = averages were carried out for an 
ensemble of 100,000 configurations, whereas for Nf — 1, 2 the average was over 20,000 
configurations. In the latter case, configurations are correlated due to the Monte-Carlo 
evolution, so that, effectively, the size of the ensemble is somewhat smaller. 

As stated in the introduction, our hypothesis is that the scale of the microscopic 
fluctuations and the scale of the secular variations are independent. The two scales are 
separated via a procedure called unfolding. This procedure is not new. It has been widely 
used in the analysis of spectra of quantum systems in terms of random matrix ensembles 

|J. It is always the unfolded spectrum that has universal properties. In the present 
case, we first obtain the average spectral density p, by fitting a smooth curve to the 
histogram of the eigenvalues of the Dirac operator. In our case we use a combination of a 
gaussian and an exponential function. This function is then used to generate the unfolded 
spectrum {A^} from the original spectrum {\ n } according to 

\ u n = [ Xn p(\)d\. (14) 
Jo 

The average spacing of the unfolded spectrum is unity. 

In Figs. 1 and 2 the solid lines show histograms of the unfolded spectrum {A^} 
for N c — 2 and N c = 3, respectively. These results are compared to the microscopic 
spectral density for the chGUE (dashed curve) and the chGOE (dotted curve). Results 
are given for 0, 1 and 2 flavors. In all cases, we find that the spectrum for N c = 2 is 
described by the chGOE, whereas the spectrum for N c = 3 is described by the chGUE. 
The only discrepancy is that the oscillations for (5 = 2 are not reproduced several level 
spacings away from zero. Since the positions of the peaks give the average location of 
the eigenvalues, this implies that level fluctuations are larger than expected according to 
random matrix theory. At this moment it is not clear whether this effect is a finite size 
effect that will go away in the thermodynamic limit. For more than two flavors the gap 
near zero becomes so big that the microscopic and the macroscopic scales can no longer 
be separated unambiguously. 

6. Discussion 
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We have calculated the spectrum of the Dirac operator for an instanton liquid ap- 
proximation to the QCD partition function for both two and three colors. For 0, 1 and 
2 flavors, we have have found that the spectrum near zero virtuality is different in both 
cases. In all cases the microscopic limit of the spectrum, obtained in the thermodynamic 
limit by rescaling the eigenvalues ~ V4, can be described by chiral random matrix theory 
with only the symmetries of the Dirac operator as input. For two colors the Dirac opera- 
tor is real and the appropriate ensemble, the chiral orthogonal ensemble, has real matrix 
elements. For three colors the Dirac operator is complex. Now the appropriate ensemble, 
the chiral unitary ensemble, has complex matrix elements. This confirms our universality 
hypothesis that the spectrum near zero virtuality is entirely determined by symmetries. 
In a sense, this is natural because this region of the spectrum carries information on the 
very long wavelength excitations which are determined by the global symmetries of the 
QCD partition function. Indeed, the microscopic correlation functions are consistent with 
the static limit of the QCD partition function. 

As has been argued in [|12|], there is a third universality class, where the matrix elements 
of the chiral random matrix ensemble are quaternion real. This case is realized for fermions 
in the adjoint representation. It would of great interest to subject this case to numerical 
studies as well. 
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Figure Captions 

Fig. 1. The microscopic spectral density ps(z) versus the virtuality z measured in units 
of the average level spacing. The results are for the color gauge group SU(2). The full line 
shows the histogram of our numerical calculations. Analytical results for the chGOE and 
the chGUE are represented by dotted and dashed curves, respectively. For the number of 
flavors we refer to the labels in the figure. 

Fig. 2. The microscopic spectral density ps(z) versus the virtuality z measured in 
units of the average level spacing. The results are for the color gauge group SU(3). Further 
explanation can be found in the caption of Fig. 1. 
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